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SUMMARY 

This report presents a preliminary analysis of the mathematical modelling of eddy 
current effects in a large-gap magnetic suspension system. It is shown that eddy currents can 
significantly affect the dynamic behaviour and control of these systems, but are amenable to 
measurement and modelling. A theoretical framework is presented, together with a comparison 
of computed and experimental data related to the Large Angle Magnetic Suspension Test 
Fixture at NASA Langley Research Center. 


1.0 INTRODUCTION 

In order to explore and develop technology required for the magnetic suspension of 
objects over large ranges of orientation, a small-scale laboratory development system, the 
Large Angle Magnetic Suspension Test Fixture (LAMSTF) has been constructed at NASA 
Langley Research Center. Possible applications for magnetic suspension systems of this general 
class include space payload pointing and manipulation, microgravity vibration isolation and 
wind tunnel model suspension [1], An important objective of this particular project is to 
investigate the dynamic modelling of large-gap magnetic suspension systems, so that future 
systems can be designed with higher confidence levels. 


2.0 HARDWARE DESCRIPTION 

The general configuration [2] is illustrated in Figures 1,2. An array of five, room 
temperature, copper electromagnets are equally spaced on a 13.77 cm radius. The coils are 
wound with 509 turns of AWG 1 0 enamelled copper wire on bakelite forms, with mild steel 
cores. The electromagnets are mounted on an aluminum plate 1.27 cm thick. Each 
electromagnet is driven by a transistor switching power amplifier, rated at ± 150 V and ± 30 
A continuous, with four-quadrant operation. 

The suspended element consists of 16 wafers of Neodymium-Iron-Boron permanent 
magnet material, each 0.851 cm in diameter and 0.3175 cm thick, epoxied into an aluminum 
tube, 5.32 cm long and 0.9525 cm outside diameter. The total mass of the suspended element 
is 22.5 grams and the moment of inertia about tranverse axes is 5.5 x 10" 6 kg.m 2 . The 
direction of magnetization is along the axis of the cylinder, which is horizontal when suspended. 
The nominal magnetization is 954,930 A/m (1.2 Tesla), although measurements have indicated 
a slightly lower working value. The suspension height is 0.1m, measured from the axis of the 
suspended element to the top plane of the electromagnet conductor. 



The position sensing system consists of multiple light beams, arranged in the vertical 
and horizontal planes, partially interrupted by the suspended element. The light sources are 
miniature infra-red light-emitting diodes, intended for use with fiber-optics, with collimating 
lenses added. The light receivers are matching infra-red phototransistors, with focusing lenses 
added. The complete sensor system is mounted on a framework which can be rotated by hand 
about a vertical axis. A schematic diagram of the sensor assembly is shown in Figure 3. 

Several different control systems have been developed and demonstrated, including a 
simple analog version with phase-advance (“lead") compensation [2], a first-generation digital 
controller, generated using the bilinear (Tustin's) transformation [3,4], a decoupled PD 
controller [5], and LQR and LQG designs [4], 


3.0 MODELLING OF EDDY CURRENT EFFECTS IN LAMSTF 

3.1 Introduction 

Whenever a time-varying magnetic flux penetrates a conducting medium, eddy-currents 
will be generated. In the case of LAMSTF, the principal time variation is due to the necessary 
control forces and torques being generated by fluctuating electromagnet currents, since the 
system is open-loop unstable. In the original design, eddy-current circuits were deliberately 
introduced in three main areas, as illustrated in Figure 4 : 

1) Position sensor structure, 2) Electromagnet cores, 3) Aluminum baseplate 

This was done so that it would be necessary to measure, analyze and model the eddy 
current effects, rather than attempting to avoid their influence, as is the usual practice. The fact 
that stable suspension was initially achieved rather easily [2] was taken to indicate that the eddy 
current effects were not very significant. However, a consistent discrepancy was discovered in 
the dynamic behaviour in the “pitch" degree-of-ffeedom, illustrated in Figure 5 [3,4], In 
consequence, an analysis and modelling effort was undertaken. 


3.2 Simplified Eddy Current Modelling 

A simplified analysis can be employed to assess the effects of eddy currents in 
LAMSTF. Two important initial assumptions are made: 

i) Some a priori knowledge of the geometry of the eddy current circuit 

ii) The circuit geometry is independent of frequency 

The first assumption might require that the eddy currents be constrained to flow around 
well-defined paths, such as the position sensor structure, rather than through large plates or 
shells of conducting material. Alternatively, the circuit geometry must be relatively simple and 
predictable. The second assumption requires that the “skin depth" be much greater than the 
local material thickness. The skin depth is given by the following formula [6,7 etc.] : 


6 = v ( — 2 —) 

\ Vo Mr O U> ) 


or 


V 


( 


J 


-(l) 



- where 5 = Skin depth (m), p = Permeability (Hm' 1 ), p= 1 la = Resistivity (0), 
<jj = angular frequency (s' 1 ). In the case of LAMSTF, the natural frequencies of the suspended 
element are rather low, of the order of 10Hz or less. For an aluminum conductor, the value at 
10Hz would be around 28mm, much greater than the typical material thickness in LAMSTF. 
The only exceptions are the iron electromagnet cores, which will be discussed at some length 
later. 

If both of the above assumptions are satisfied, the resulting model corresponds to that 
occasionally described in literature as the Single Time Constant Model. The derivation 
resembles the analysis of a transformer with a shorted secondary, as illustrated in Figure 6 : 

V = IR + Lfl+L m ,2u+L m; ^ +... -(2) 
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- where R^, are the resistance and inductance of the n'th eddy current circuit and L mn is the 
mutual inductance between the primary (the electromagnet coil) and the eddy current circuit. 
Note that mutual inductances between multiple eddy current circuits are neglected. The 
terminal characteristics of the primary (driven coil) can be found by combining these two 
equations : 
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One special case is of interest here. Suppose that : 
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but if Rc,-» 0 or s -* oo : 
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This confirms that a non-dissipative (reactive) secondary effectively “turns off' part of 
the primary inductor. Continuing, the field components generated (at the suspended object) 
can be expressed as : 


B i =K i I + K,,U,+- = Kjl ( 1 - i^X, --) 
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- where Kj, K^ are constants representing the j'th field component generated at the 
suspension location by the electromagnet and the n'th eddy current respectively. Now the factor 

will, in general, be different for each field component, that is each individual eddy current 
will affect each field component by a different proportion, compared to the main electromagnet. 
Therefore the eddy current effects in a system involving several electromagnets and eddy 
current circuits should be represented as follows : 

[Bj] = [KJ [I] + [KjJ [Ie] - (8) 

- where [Bj] = (B x B y B z Bxx ••• ) T , [I]=(Ij I 2 ••• ) T , [Kj] is a rectangular matrix of 
field coefficients and [KjJ is similar, though possibly of differing dimension. It is presumed that 
[Ie] can be derived from [I], following equations 2,3. 

Alternatively, if the eddy current circuit has similar geometry to the primary (for 
example the induced current in electromagnet cores), it can be argued that the relative effect on 
all field and field gradient components at the suspended object will be similar. In this case, the 
representation can be considerably simplified by invoking a false current as shown : 

I'"(l - ie, (*%,%) - - ) 1 . where Bj = Kjl' -(9) 

It should be noted that the change in electromagnet terminal characteristics and the 
change in field at the model are two separate effects and should be modelled as such. 

3.3 Alternative Representations - The Circuit Model 

The preceeding analysis was cast in simple terms, based on the concept of flux linkages 
between electrically separate circuits. This corresponds to the physical situation, and is a 
reasonably comprehensible representation for non-specialists. A significant amount of previous 
literature does exist in the area of eddy current modelling, but is mostly cast in terms more 
familiar to electrical engineers, namely circuit impedances [8,9,10 etc.]. So as to clarify the 
equivalences between different approaches, an alternative representation will now be 
established. 

Considering a single coil, wound on a dissipative core, illustrated in Figure 7, the 
equation governing the magnetic excitation is as follows : 

MMF = G e f t + i <f> 
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- where MMF is the magnetomotive force (ampere-turns) driving the magnetic circuit, 
G e is the electrical conductance (inverse of resistance) and P is the magnetic permeance. The 
magnetic impedance is defined as the complex ratio of MMF to flux rate : 
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Now using equation 10 and assuming sinusoidal excitation as usual : 
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Now the concept of the “gyrator" is introduced, wherein the magnetic impedance can be related 
to the electrical impedance of the exciting coil : 



The terminal characteristics of the coil can now be established : 

i _ j_ _ i _ i 
V Z T R+Z« r+J±£- 


-(13) 


-(14) 


This is equivalent to equation 4 above, as can be verified with the following substitutions and 
manipulations : 

L e , = ^2 (assume eddy current geometry resembles driving coil) -(15a) 

Lm, =^(LL e| ) = | (perfect flux linkage to eddy current circuit) - (15b) 

Rc, = (definition) -(15c) 

L = N 2 P (definition) - (1 5d) 

Note that equations 15a and 15d indicate that P = L e , in this representation. 


4.0 EXPERIMENTAL MEASUREMENTS 
4 . 1 Determination of Parameters 

The question now is, can the parameters K^, L^, R^ and L mn be estimated and/or 
measured with sufficient accuracy? First the problem of estimation is addressed. 

Calculations have been carried out using the finite element computer code VF/GFUN, 
by Vector Fields Inc.. It should be noted that this code is magnetostatic and has no capability 
for direct eddy current calculations, although such codes are available (for instance ELEKTRA, 



by the same supplier). Instead, the code is used to calculate flux linkages, hence inductances, 
using : 
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VF/GFUN calculates the field on a grid representing the linkage plane of the eddy 
current circuit. The field normal to the plane is then numerically integrated (by the OPERA 
pre- and post-processor) to yield the flux linkage term. Figure 8 illustrates the general 
arrangement. The calculation of the Ke jn terms is straightforward. 

By way of example, a series of calculations has been made for a single LAMSTF 
electromagnet with a representation of one part of the position sensor assembly mounted on the 
same axis, as illustrated in Figure 9. The required parameters were predicted (or previously 
measured) to be : 

L = 0.0275 H R = 0.74 ft 

L e = 6.69 x 10' 7 H Re = 2.243 x 10 4 Q 

L m = 1.0998 x 10' 5 H 

K z = 3.495 x 10 4 T = 4.369 x 10 -6 T 

Incorporating these values in equation 7, and examining the axial (z-axis) field 
component, gives : 


B z = K Z I 
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It is seen that the reasonant frequency of this eddy current circuit is around 53Hz, 
significantly higher than LAMSTF open-loop natural frequencies, but still well within the range 
of interest. 


4.2 Experimental Verification 

Actual measurements of the current to field transfer function, corresponding to equation 
7, were made with an experimental set-up as described above, and later with LAMSTF. Field 
components were measured with a F.W. Bell Model 9903 Hall-effect gaussmeter. 
Electromagnet currents were measured using a current shunt. The transfer function was 
measured directly with a Schlumberger Model SI 1250 analyzer, with sine-sweep excitation. 
The results for an air-cored electromagnet with no eddy current circuits are shown in Figure 
10, and are taken to represent the probe + instrument + data acquisition system response. 
These results are subtracted from all subsequent measurements. Figure 1 1 shows the measured 
transfer function for B z , together with the predictions from equation 17. The agreement is 
thought to be satisfactory. The values of most parameters can be adjusted (refined) to give a 
better agreement, as shown in Figure 12. The only significant residual discrepancies are seen to 
occur at higher frequencies where the validity of the Single Time Constant Model is 
questionable. 



The adjusted form of equation 1 7 is : 


b, = kj ( i - ) -07b) 

4.2. 1 Parameter Adjustment Procedure 

If equation 1 7 is written : 

Bj = K jl ( l - life ) - 08 ) 

- then it is easily shown that the maximum phase lag occurs at a frequency given by : 
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- that the maximum phase lag is : 

‘T-’max = Tan ^ 27cb^abj ) " @0) 

- and that the infinite frequency gain attenuation is : 

dtw> = 20 logio (l-§) -(21) 

In most cases, the maximum phase lag and the frequency for maximum phase lag are the easiest 
to apply, and are sufficient to specify the parameters a and b, using : 

an(J b=f+v(^ + ^-) -(22) 

4.3 More Complex Cases 

If the electromagnet is mounted on the aluminum plate, a second eddy current circuit is 
added; when the iron core is inserted, a third is added. Figures 13-15 show the comparison 
between experimental and computed responses. Again, the agreement is fair, although capable 
of improvement by refinement of parameter estimates, also shown in Figures 13-15. Note that, 
even if refinement of parameters is undertaken, the model does not correctly predict the high 
frequency behaviour, particularly where the iron core is present. This is due to the iron core 
becoming highly dissipative at these frequencies, due to its small skin depth. 

An additional series of calculations and measurements has been made for the vertical 
field component generated at the centroid of the suspended element due to a LAMSTF 
electromagnet at the design location, i.e. with the sensor ring off-axis relative to the 
electromagnet. The geometry is illustrated in Figure 16 and typical measured responses is 
shown in Figure 17. It is clear that the model is progressively less suitable for progressively 
more complex geometries. 



Measurements have also been made with the full LAMSTF sensor frame in place. 
Examples are shown in Figure 18. Due to the geometrical complexity, with multiple 
interlocking eddy current paths, it has not yet been possible to compute resonable estimates by 
the methods shown. It should be noted, however, that another assumption inherent in the form 
of the model previously chosen has been violated, that is the lack of interaction between 
separate eddy circuits. If these circuits are physically and electrically connected, this is clearly 
not reasonable. 


4.4 Terminal Characteristics 

It appears to be possible to experimentally estimate certain important parameters 
without direct measurement of magnetic fields. Figure 19 shows a comparison of measured 
and computed terminal characteristics for the single LAMSTF electromagnet mentioned above. 
The agreement is not perfect, but sufficient to validate the approach and can, of course, be 
improved by adjustment of parameters. 


5.0 ADVANCED EDDY CURRENT MODELLING 


5.1 Introduction 

The discrepancies between the modelled and measured responses with the iron core 
present are confined to the higher frequencies. They are manifested as reduced intensity and 
increased phase lag of the magnetic field, compared to the simple model described above. The 
simple analytic model cannot be adjusted or manipulated so as to satisfactorily model this 
behaviour. However, by studying some more simplified geometries, where analytic solutions 
are available, insight into the correct form of a more elaborate and realistic model can be 
gained. The cases to be studied are the magnetization of an infinite sheet by an in-plane 
uniform applied field and the magnetization of an infinite cylindrical bar by a uniform axial field. 


5.2 Magnetization of an Infinite Flat Sheet 

This is a classical problem, whose solution can be found in numerous texts [6,7,11,12] 
and which is illustrated in Figure 20. The outline is shown here for reference. The governing 
equation is the 1 -dimensional diffusion equation : 
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This can also be written is terms or H z , or, in an alternative form, in terms of J x . The general 
solution to this equation is : 


B z = K,e Qy + K 2 e“ Qy 
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From the boundary conditions of B = B s @ y = ± b, we find : 


Ki = K 2 


e ob_j_ e -ob 


So that : 



This is plotted as Figure 21. Integrating the total flux acrosss the conductor : 
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Now expanding the hyperbolic functions as infinite series : 
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Where : 

(ab) 2 = jw (a/i 0 /i r b 2 ) -(30) 

It is recognized that equation 29 represents a pole-zero sequence, with a d.c. response of 2B s b. 
The obvious approach to modelling is to approximate using the first few terms of the series. 
Figure 22 shows the response for various numbers of terms. It is seen that large fluctuations in 
response occur at higher frequencies with larger numbers of terms. This is due to the 
generation of lightly damped, closely spaced poles, coupled with series truncation error. In 
practice, only a few terms would be employed in any case, to avoid excessive complexity in the 
model. By inspection, a four- or five-term series seems to be a reasonable choice. 

It should also be noted that the limiting rate of roll-off (i.e. at high frequencies) is 
lOdb/decade and that the limiting phase lag is 45°. Both of these values is exactly half the 
corresponding value for a simple first-order pole. 

5.3 Magnetization of an Infinite Cylindrical Bar 

The magnetization of an infinitely long cylindrical bar in a uniform, axially applied field 
is again a classical problem whose solution can be found in a number of texts [6,7 etc.]. The 
geometry is shown in Figure 23. The governing equation is the diffusion equation : 

V 2 B = Wr f -(31) 

In cylindrical coordinates, for an axisymmetric problem, this reduces to : 



For the case of sinusoidal excitation, we assume B x = B(r) e’ 1 ** and the spatial part of the 
problem reduces to the form of Bessel's equation : 

= » -(33) 

The solution of this equation is of the form : 

B = Aj J 0 (ar) + A 2 Lo(ar) - (34) 

- where J 0 and are zero-order Bessel functions of the first and second kind respectively, and 
a has the same definition as previously. Since a finite field is required on the axis (r = 0), the 
coefficient A 2 must be 0. With the boundary condition of B = B s at r = a, the solution 
becomes : 
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- which can be expressed as an infinite series : 
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This is plotted in Figure 24. Again integrating the total flux across the conductor : 
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This is shown in Figure 25. 


6.0 SCALING of EDDY CURRENTS 

Whenever sub-scale tests are contemplated, it is necessary to understand how the 
phenomena under investigation scales with linear dimension. For instance, in this case, sub- 
scale tests may well be required to validate predicted eddy current models, or even to provide 
design values for some parameters. 

For the case of conducting ring of given geometry and scale S, as shown in Figure 26, it 
is immediately apparent that : 

Rap| a ^ Loc N 2 S Therefore - ^ oc S 2 - (38) 

Of course the proposed dynamic models (for instance see equation 4) are composed of terms of 
the form (R + ja;L). Now recalling the skin depth formula : 



In either case it is clear that dynamic similarity can be achieved if the excitation frequencies are 
scaled as follows : 
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7.0 DISCUSSION 

The eddy current models proposed appear to be satisfactory in the case of large eddy 
current circuits in conducting, magnetic, or non-magnetic material. Relatively simple 
computations are capable of providing reasonable estimates of important parameters, with the 
option of refinement based on measurements of magnetic field or electromagnet terminal 
characteristics. The most difficult case is seen to be that of eddy currents in magnetic material, 
due to the lower skin depths. However, discrepancies at frequencies well outside (above) the 
controller bandwidth would be of no particular consequence in this case. 

The simple models are not adequate to fully describe the eddy currents in the LAMSTF 
position sensor structure, due to its geometric complexity. This leads to difficulties in 
parameter estimation and strong coupling between multiple eddy current circuits. 


8.0 CONCLUSIONS 

Simple models for the effect of eddy currents in the metallic structure of LAMSTF has 
been proposed and validated by experiment. However, it has not yet been possible to fully 
describe the eddy currents in the position sensor framework due to the geometric complexity 
involved. 

Eddy currents have been shown to seriously affect field and field gradient components 
in the frequency range of interest, such that they must be incorporated into a system dynamic 
model if modem control synthesis techniques are to be fully successful. 
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Figure 3 - Schematic Diagram of Optical System 
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Figure 4 - Schematic Diagram of Eddy Current Circuits in LAMSTF 
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Figure 5 - Initial Discrepancy in Dynamic Response in Pitch (0 y ) 



Figure 6 - Schematic of Circuit Model for Primary to Secondary Coupling 





Figure 8 - Illustration of OPERA Flux Linkage Computation (B z shown) 




Figure 9 - Simplified Eddy Current Test Set-Up 
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Figure 10 - Gaussmeter Probe Response 
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Figure 1 1 - Axial Field for Air-Cored Electromagnet and Dummy Sensor Ring 
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Figure 12 - Axial Field with Adjusted Parameters 
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Figure 1 3 - Axial Field with Alloy Plate Added 
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Figure 14 - Axial Field with Iron Core Added 






Figure 1 5 - Axial Field with Iron Core and Alloy Plate Added 






Figure 16 - Off-Axis Test Geometry 
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Figure 17 - Vertical Field with Off-Axis Test Set-Up 
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Figure 19 - Electromagnet Terminal Characteristics 




Figure 21 - Field Penetration for Infinite Flat Sheet Problem 




Figure 22 - Integrated Flux for Infinite Flat Sheet Problem 
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Field Penetration into Cylindrical Bar 



Figure 24 - Field Penetration for Infinite Cylindrical Bar Problem 
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Figure 25 - Integrated Flux for Infinite Cylindrical Bar Problem 







